DOUBLE SCALING LIMIT IN RANDOM MATRIX MODELS AND A NONLINEAR 
HIERARCHY OF DIFFERENTIAL EQUATIONS 
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Abstract. We derive the double scaling limit of eigenvalue correlations in the random matrix model at 
critical points and we relate it to a nonlinear hierarchy of ordinary differential equations. 



1. Introduction 
We consider the unitary ensemble of random matrices, 

dfi N {M) = Z N X exp (— iVTr V{M)) dM, Z N = 



exp (-N Tr V(M)) dM, 



(1.1) 



on the space Hn of Hermitian N x N matrices M — (My) 1< . . <N , where V(x) is a polynomial, V(x) = 
v p x p + Vp^ix?' 1 + . . ., of an even degree p w ith v p > 0. The ensemble of eigenvalues A = {Aj, j = 1, . . . , N} 



of M is given then by the formula (see e.g. [Mehj, JTW |) 



dfj, N (X) = Z- 1 exp (-NH N (X)) dX, Z N = exp (-NH N (X)) dX, 

Ja n 



(1.2) 



where An is the symmetrized 



pJV 



A 



bN 



/S(N), and 



A? 



H N (X) = -- ]T l°g|Ai-A fc |+5^V(A,-). (1.3) 

i<i<fc<w j=i 

Let dv^ix) — pN{x)dx be the distribution of the eigenvalues on the line, so that for any test function 
<p{x) € C °°, 

N 



djjL N {X) 



<p(x) dv^{x) 



As JV — > oo, there exists a weak limit of dv^{x)^ 



di/ 00 (x) = lim di>N(x) . 

N—>oo 



(1.4) 



(1.5) 



To determine the limit (cf. [BIPZ|, [DGZ], and others), consider the energy functional on the space of 

(1.6) 



probability measures on the line, 

I{dv(x)) = - // log |a: - y\dv{x) dv{y) + / V{x)dv(y) 



Then Hn(X) in fll.q ) can be written as 

H N (X) = NI(di/(x; A)) , 
where dv{x\ A) is a discrete probability measure with atoms at Aj, 

1 W 

dv(x; X) — — 5 (a; — Aj) da; . 



(1.7) 



(1.8) 



j'=i 
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Hence, 

dfi N (X) = Z N X exp (-N 2 I(dis{x; A))) dX . (1.9) 

Because of the factor N 2 in the exponent, one can expect that as N — *■ oo, the measures g?^at(A) are localized 
in a shrinking vicinity of an equilibrium measure dv eq (x), which minimizes the functional I{dv(x)), and 
therefore, one expects the limit (1.5) to exist with dvoo(z) = dv oq (x). A rigorous proof of the existence and 
uniqueness of the equilibrium measure, its properties, and the existence of limit (1.5) with dvoo(x) — dv cq {x), 
was given in [BPS] and [Joh|. 

The equilibrium measure dv eq (x) is supported by a finite number of segments [aj, bj], j — 1, . . . , q, and 
it is absolutely continuous with respect to the Lebesgue measure, dv eq (x) = p(x)dx, with a density function 
p(x) of the form 

i i 
Pi*) = ^-Hx)R 1 l 2 (x), R{x) = H(x aj )(x - bj), (1.10) 

where h(x) is a polynomial of the degree, deg h = p — q — 1, and R 1 J_ 2 {x) means the value on the upper cut of 
the principal sheet of the function R 1 / 2 (z) with c uts on J. The equilibrium measure is uniquely determined 
by the Euler-Lagrange conditions (see |DKMVZ |): for some real constant I, 



log \x — s\ dvcq(s) — V(x) — I, for x € U? =1 [oj, bj] , 
log \x — s\ dv eq (s) — V(x) < I, for Uj =1 [aj,bj 



Equations ( 1.10 ), (|l.ll| ) imply that 



uj(z) 



V'(z) h(z)R 1 / 2 (z) 



where 



L)(z) 



In addition, (1.11) implies that 



p(x) dx 



h{x)R l ' 2 {x) 



z^ + Oiz- 2 ), z 



oo . 



dx = 0, j = 1,. .. ,q- 1, 



which shows that h{x) has at least one zero on each interval bj < x < ; j = 1, . 
we obtain that 



V'(z) = Pol h{z)R 1/2 {z) 



Res 

2 — OO 



h{z)R 1/2 {z) 



and 



h{z) 



Pol 



V'(z) 



(1.11) 
(1.12) 

(1.13) 
(1.14) 

(1.15) 

,q- 1. From ( |1.13| ) 
(1.16) 

(1.17) 



.i?V2( z )_ 

where Pol[/(z)] is the polynomial part of f(z) at z = oo. The latter equation express es fe(z ) in terms of 
V(z) and the end-points, <zi, , q q , 6 q . The end-points can be further found from ( 1.16| ), which gives 

q + 1 equation on ai, . . . , b qi and from (1.15), which gives the remaining q — 1 equation. 
The equilibrium measure <iz/ eq (:c) is called regular (otherwise singular), sec [DKMVZ], if 

h(x)^0 for xe U| =1 [aj,6j] (1.18) 

and 



2 / log |x — s| di/ oq (s) — ^(x) < I, for sc G 



(1.19) 



The polynomial F(x) is called critical if the corresponding equilibrium mesaure dv eq (x) is singular. To 
study the critical behavior in a vicinity of a critical polynomial V(x), one embeds V(x) into a parametric 
family V(x; t), t = {t\, . . . , t r ), so that for some t c , V(x; t c ) — V(x), and the problem is then to evaluate the 
asymptotics of eigenvalue correlation functions as t — > t c . The number of parameters r depends, in general, 
on the degree of degeneracy of the equilibrium measure <i^eq(a;). 
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In this paper we concern with the critical behavior for the polynomial V(x) such that the corresponding 
equilibrium measure is supported by the segment [—2, 2], with a density function of the form 



p(x)^Z- 1 (x-c) 2m y/A- 



(x- c) 2m ^A-x 2 dx, 



(1.20) 



where — 2 < c < 2 and m = 1,2,.... The choice of the support segment is obviously not important, because 
by a shift and a dilation one can reduce any segment to [—2, 2]. The parameter m d eterm ines the degree of 
degeneracy of the equilibrium measure at x — c. Observe that when c ^ 0, density ( 1.20 ) is not symmetric. 
Our results are summarized as follows. We are interested in two problems: 

(1) The singularity of the infinite volume free energy at the critical point. 

(2) The double scaling limit, i.e. the limit of rescaled correlation functions as simultaneously the volume 
goes to infinity and the parameter t goes to t c , with an appropriate relation between t — t c and the 
volume. 

Case m = 1. Free energy. We evaluate the derivatives in T of the (infinite volume) free energy F(T), 
where T > is the temperature, and we show that F (T) can be analytically continued through the critical 
value T c both from below and from above of T c . In addition, we show that F(T) and its first two derivatives 
are continuous at T = T c , while the third derivative has a jump. This proves that at T = T c the phase 
transition is of the third order. It gives an extension of the result of [GW| where the third order phase 
transition was shown for the case of a symmetric critical V{x) in the circular ensemble of random matrices. 

Case m=l. Double scaling limit. The key problem here is to derive a uniform asymptotic formula for 
the recurrence coefficients of the corresponding orthogonal polynomials. The double scaling limit descibes a 



transition from a fixed point behavior of the recurrence coefficients to a quasiperiodic behavior (cf. [DKMVZ ] 
and [ BDE ] ) , and the problem is to derive a uniform asymptotic formula for the recurrence coefficients in the 
transition region. We show that under a proper substitution, the recurrence coefficients are expressed, with 
a uniform error term, in terms of the Hastings-McLeod solution to the P ainle vc II d ifferential equation. In 
the symmetric case (c = 0) our solution reduces to the one obta in ed in f5|f], [jPeStl . For a rigorous proof 



of the double scaling asymptotics in the symmet ric ca s e sec | BI2|, [BDJ j. The both latter papers are based 



on the Rieman n-Hilb er t ap pr oach , developed in | FIK | , [ BI1 |, | DKMVZ]] . It is worth mentioning also earlier 
physical works BKa |, [ GM |, D£ | which concern with the double scaling limit of the Painleve I type. 

General case, m > 1. Double scaling limit. We derive a hierarchy of nonlinear ordinary differential 
equations which give, under a proper substitution, the double scaling limit of the recurrence coefficients for 
all m > 1. The hierarchy admits a Lax pair of linear differential equations and it can be constructed in the 
framework of the general theory of isomonodromic deformations [IN]. Our particular hierarchy is known as 
the P ainlevc II hierarchy |Kit| and it is related to selfsimilar solutions of the mKdV equation | PeS (see also 
] Nfoo| ). 



2. Critical Behavior for a Nonsymmetric Quartic Polynomial 
Let us consider the critical quartic polynomial V c (x) such that 



Vo(x) = ^(x 3 -4 Cl x 2 + 2c 2 x + 8 Cl ), T c = l + 4c 2 ; V c (0) = , 



where we denote 



Ck = cos kite , Sfc = sin kite. 



This corresponds to the critical density 



{x) = ^ F {x-2c 1 ) 2 ^x~ 2 - 



(2.1) 
(2.2) 
(2.3) 



Observe that < e < 1 is a parameter of the problem which determines the location of the critical point, 

-2 < 2ci = 2cos?re < 2. (2.4) 



Equation (1.13) reads in this case as 



V' c (z) (z - 2c 1 ) 2 V / 7 r ^4 



2T„ 



(2.5) 
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The correlations between eigenvalues in the matrix model are expressed in terms of orthogonal polynomials 
Pn{x) = x n + ... on the line with respect to the weight e~ NVc ^ (see e.g. |Mch], | TW |). Let 

4, n ( x ) = -^P n ( x )e- NV ^/ 2 , n = 0,l,..., (2.6) 

be the corresponding psi- functions, which form an orthonormal basis in L 2 . They satisfy the basic recurrence 
relation (see e.g. [fBzef ), 



h„ 



(2.7) 



and the differential equation, 
1 V'(x) n 1 1 1 

(2.8) 

(2.9) 
(2-10) 



The compatibility condition of equations (2/7), (2^) leads to the string equations, 

T ^ T \ = ll + 7„-l + 7n+l + Pi + $1-1 + PnPn-l ~ 4c X (/3„ + /3„_i) + 2c 2 , 



■N % 



= V c '(/?„) + 7 ^(2/3„ + /3„_i) + 7 " +1 (2/3n + /3„+i) - 4 Cl ( 7 * + 7r 2 l+1 ) 



To study the critical asymptotics we embed V c {x) into a parametric family of polynomials. To that end 
for any T > we define the polynomial 



V(x;T) = ±V(x), 



(2-11) 



where V(x) is such that 

V'[x) = x 3 - A Cl x 2 + 2c 2 x + 8 Cl , \/(0) =0. (2.12) 

Then V^'(ar) = V^'(j:;r c ). We call T temperature and T c critical temperature. Denote AT = T — T c . Let 
p(x;T) be the equilibrium density for the polynomial V(x;T). Equation ( 1.13| ) reads in this case, 

f p{x;T)dx _V'(z) h{z;T)R 1 / 2 {z;T) 
h{T) 

where h(z; T) is a monic polynomial in z. 



,(z;T) 



z — x 



2T 



2T 



(2-13) 



Free energy near the critical point. The (infinite volume) free energy is defined as 



F(T) 



-T lim —\nZ N (T), Z N {T) 

N— >oo iv 



exp 



N 



Tr V(M) dM 



(2.14) 



We will show that at T = T c , F(T) is not analytic. To evaluate the type of nonanalyticity at T = T c , 
consider the function 



2 d fF(T)\ 



It can be evaluated as 



= lim — 



%— [ — Tr V(M) exp (-— Tr V(M)) dM. 
(T) Jn N N \ T V 7 



(2.15) 



Fx{T)= [ 

Jj(T) 



V(x)p(x; T)dx 



1 



AniT 



V(z)h{z;T)y/R(z;T) dz 



/ V(z)w(z;T)dz, 
2m J c 



(2.16) 



where C is any contour with positive orientation around J(T), the support of equilibrium measure. Observe 
that the limits in (2.14), (|2.15 ) exist for any polynomial V(x), due to the weak convergence of eigenvalue 
correlation functions (cf. |Joh| ). On the contrary, for the second derivative of F(T; N) = —(T/N 2 ) In Z^(T) 
in T, the convergence F"(T;N) — > F"(T) does not hold if the equilibrium measure has two cuts or more, 
because of quasiperiodic oscillations of F"(T; N) as a function of N (see [BDE ). From (2.16) it follows that 
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since uj(z;T) is continuous on C in T at T = T c , F\(T) is continuous as well. Therefore, F'(T) is continuous 
at T = T c . Consider F"{T). 



Second derivative of the free energy. From ( 2.16 ), 

d[TFx{T)] 



dT 



(2.17) 



For T > T c , the equilibrium measure corresponding to V(x : T) is supported by one cut [a,b] and the 
equilibrium density is written as 



P(x; T) = -L [(a- _ C )2 + d 2 ] V(b-a:)(a;-a) , 
where a = a(T), b — b(T), c = c(T), <i = d(T). In the one-cut case we have the equation, 



see Appendix below, hence 



This implies that 



djTF^T)} 
dT 



d[TFi(T)} 



V(z) 



2ni Jc y/( z -a)(z- b) 



: dz. 



dT 



1 



T=T, 



+ 2ni J c yf 



V{z) 



: dz. 



We should mention here another useful formula valid for T > T c (see [DGZ|, |Eyn|): 

d 2 [TF{T)] _^ 2ln b-a 



(2.18) 
(2.19) 
(2.20) 
(2.21) 

(2.22) 



dT 2 " _ 4 

For T < T c , the equilibrium measure corresponding to V(x) is supported by two cuts [ax, bx] and [aa, 62] 
The equilibrium density is written in this case as 



1 



P( x l T ) = 7T^ ( x ~ c ) V ( & i - x )i x ~ ai){b 2 ~x)(x- o 2 ) , 

Z7TJ 

where where ai, b\, a 2l b 2 , c depend on T and b\ < c < a 2 - In the two-cut case we have the equation, 



4f[Tu(z',T)] = Z X ° = 

dT 1 v/(z - - - a 2 )(^ - fe 2 ) 

where xq = Xq(T), b\ < xq < a 2 , is determined from the condition that 

: dx = , 



bi \J (x - ax)(x - bx){x - a 2 )(x - b 2 ) 

V(z)(z - x ) 



see Appendix below, hence 

d[TFx(T)] = J. 

dT 2ni J c [ z _ ai )(z - b 1 ){z - a 2 )(z - b 2 ) 

We have that a 2 = b\ = x n = 2c\ at T — T~ , hence 



: dz. 



d[TFx(T)} 



Thus, 



dT 
d[TFx(T)} 



dT 



T=T~ 



T=T~ 



2ni Jc y/z 2 - 4 
d[TFx(T)} 



dz. 



dT 



T=T~J 



so that F"(T) is continuous at T = T c . Consider now F"'(T). 

Third derivative of the free energy. In the one-cut case we have that 

±l( x - c ) 2 + d 2 W(b-x)(x-a) = - \ 

aJ- y(b— x)(x — a) 



(2.23) 
(2.24) 

(2.25) 

(2.26) 

(2.27) 
(2.28) 

(2.29) 
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and 



da 



db 



dT h(a)(a-b)' dT h(b)(b - a) 



h(x) = {x- cf + d 2 



see Appendix below. From (1.17) we find that 

a + b 



c = 2ci — 



d 2 = ^ (a + bf - ci(a + b) - i ab - 2 , 



(2.30) 
(2.31) 

%)U-2, .=2 = 4(d + l) 2 , h(b)\ a= _ 2i b=2 = 4( Cl - l) 2 . (2.32) 

Therefore, a(T) and o(T) are analytic at T = T+, as a solution of system ( [2.30 ) with analytic coefficients. 
Equation fl2.20| ) implies that Ti(T), and hence F(T), are analytic at T = T+. From ( |2.30D and ( |2.3l| ) we 
obtain 



and then that 



da 



T=T,7 



4(1 + Cl f 



db 
dT 



1 



T=T< 



+ 4(1- Cl )2 

The analyticity at T = T~ is more difficult . In the two-cut case we have that 



dT 



(z - c) \/(z- ai)(z - bi)(z - a 2 )(z - b 2 ) 



2(z - x ) 



where b\ < xq < a 2 solves equation ( 2.25 ), and 
da\ 4(ai — x ) 



dbo 



\/(z- ai)(z - b 1 )(z - a 2 )(z - b 2 ) 

4(6 2 - X ) 



dT (ai - c)(ai - &i)(ai - a 2 )(ai - b 2 ) ' dT (6 2 - c)(& 2 - ai)(&2 - &i)(&2 - 02) 
see Appendix below. At T = T~ this gives that 



da 1 
~dT 



T=t: 



4(1 + Cl ) 2 ' 



d&2 
dT 



4(1 -ci) 2 



Define d and S such that 



c — d + 5, a 2 = c + d + 5 . 



(2.33) 

(2.34) 

(2.35) 

(2.36) 
(2.37) 



Then d, 5 -> as T -> T" and from (|l,15|) , 



= 



(a; — c) ^(c + d + (5 — x) (2; — c + d — (5) y/ (x — a\){b 2 — x) dx 



(u + S) V ' d 2 — u 2 y/ '(u + c — a\ + S)(b 2 — c — 5 — u) du 



d 3 J (^v + \fl - v 2 y/(vd + c - ai + 5){b 2 -c-S-vd)dv 



d 3 y/(c- ai )(b 2 - c) 



l ( [ai + b ^ 2c \ f v 2 VT^d v+ s - f VT^ dv 

2 (c- oi)(6 2 - c) J_ x d J_! 



0[d i + — 



d 3 y/{c- ai)(b 2 - c) 



it (ai + 6 2 - 2c) 
16 (c- ai)(6 2 - c) " ' 2 d 



r/F~~ + O ! -A' ! - 



hence 



1 ^+b2-2c) d2 + oid4) 



8 (c- ai )(6 2 -c) 

By the implicit function theorem, S is an analytic function of d 2 at d = 0. Equation ( 1.17 ) gives that 

Oi + b 2 1 



(2.38) 
(2.39) 



c = 2c x - 



4 



1 5 

d 2 = 2(1 + a)dai - 2(1 - ci)d6 2 - r <5 2 - - (da 2 + dfe; 2 .) - ~ daido 2 

2 8 4 

dai = ai + 2, do 2 = —2 + 6 2 . 



1 



(2.40) 
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Therefore, from fl2.36Q and fl2.39|) we obtain that as T -> T ( 



d 2 = --^AT + 0(AT 2 ), S = --^AT + 0(AT 2 ); 



c = 2ci 



3ci 
16^ 



AT + 0{AT 2 



^±^ = 2ci _i£L A r + o(AT 2 ) 



16s? 



Define now Sq such that 

61 = x — d + S , a 2 = x + d + So 
Then, similar to ( |2.38| ), we derive from ( [2.25 ) that 

: 



(2.41) 
(2.42) 



x — xq 



o-d+So \fix~o + d + So - x){x - x + d - S ) ^{x- ai)(b 2 - x) 



d.r 



y/{c-ai)(b 2 - c) 



7r (ai + 6 2 — 2c) 
2 ( c _ ai )(6 2 -c) 



d + it ^ + O [ d 



hence 



so that as T -> T~, 



2 (c- ai)(6 2 - c) 



S = ^AT + O (AT 2 ) 



(2.43) 
(2.44) 



(2.45) 



By the implicit function theorem So is an analytic function of d 2 at d = 0. Using (2.41) we obtain that 



x Q = 2 Cl -—^AT + 0(AT 2 ) 
16 sj 



By (2.37) and (2.42), 



xo = c + S — S , 



hence equations (2.35) can be written as 
dai 4(cti — c — S + So) 



db 2 



4(b 2 - c - S + So) 



(2.46) 
(2.47) 
(2.48) 



dT ( ai - c){d 2 - ( ai - c - S) 2 ){ ai - b 2 ) ' dT {b 2 -c){d 2 ~{b 2 -c-5) 2 ){b 2 -a 2 )' 

Observe that d 2 is an analytic function of a\, b 2 , and 5, So are analytic functions of d 2 . This gives the right 
hand side in (2.4S) as analytic functions of 01, b 2 and hence a\, b 2 are analytic as functions of T at T = T~ . 
Set 

bi + a 2 



Then b\ — m — d, a 2 — m + d, hence, by ( p. 26] ) 
d[TF x {T)] 1 



dT 



V{z)(z - xo) 



: dz. 



27ri J c ^(z — ai)(z 2 - 2mz + m 2 - d 2 )(z - b 2 ) 

Since Xo, m and d 2 are analytic in T at T = T~, we obtain that F(T) is analytic at T = T~ 
By (p0|), 



d^TF^T)] 



dT 2 



and by (gjSg), 

d^TF^T)} 



dT 2 



2ni 



2ni 



V(z) 



V(z) 



d 



dT V ^( Z _ a )( z _ 6 ) 



dz . 



(2.49) 
(2.50) 

(2.51) 



T=T+ 



(2 - Xq) 



dT \ y/( z - ai)(z - b 2 ) J V ' z 2 - 2mz + m 2 - d 2 



+ 



(z - x ) 



y/(z — a\)(z — b 2 ) dT \ Vz 2 — 2mz + m 2 — d 2 



dz . 



(2.52) 



T=T~ 



s 
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Observe that by (§]||), 



da 
dT 



T=T? 



dai 

It 



T=T~ 



4(1 + ci) 2 ' 



db 
dT 



T=T? 



db2 
dT 



1 



T=T' 



4(1 -d) 2 



(2.53) 



hence 



dT 2 



<P[TF X {T)] 



T=T+ 



dT 2 



T=T~ 



■5? ' 



(z - Xq) 



\Jz 2 - 4 dT \y ' z 2 - 2mz + m 2 - d 2 



dz . 



T=T~ 



From (2.41) and (2.46) we find that 



(z - X ) 



hence 



dT \ y/ z 2 — 2mz + m 2 — d 2 
d 2 [TTi(T)] d 2 [TF 1 (T)} 



dT 2 



dT 2 



T=T~ 



T=T~ 



2ni 



V(z) 



1 — C\ z + c\ 
s\{z-2 Cl ) 2 ' 

1 1 — CiZ + C 2 



(2.54) 



(2.55) 



v / i 2 ~^4 s \( z -2 Cl y 



dz 



3 + 25c 2 + 2c\ 



< 0. 



(2.56) 



Thus, F{T) is analytic both at T = T+ and T = T~, and F"'(T) has a jump at T = T c . Therefore, T = T c 
is a critical point of the third order phase transition. 

Recurrence coefficients near the critical point. The recurrence coefficients 7n , (3 n approach fixed values 
for T > T c , see e.g. ]DKMVZ|]. Namely, 



lim j n =j(T), lim /3 n =/3(T), 

n,N— >oc; >X- 



where 7 = 7 (T), /3 = /3(T) are fixed points of (|J), fl2.10|) , so that 

T-. = 3 7 2 + 3/3 2 - 8ci/3 + 2c 2 , 
= + 6 7 2 /5 - 8ci 7 2 . 

The values 7 = 7 (T), j3 — (3{T) can be expressed in terms of the end-points of the cut as 

b — a b + a 

7 = — ;— , P = 7, , 



see e.g. 



DGZ| . Therefore, by (|3|), 
1 + c 2 



7 = 1 + 



, AT + 0(AT 2 ) , /? = -— j AT + (3(AT 2 ) , AT — > H 

8s* 2sj 



(2.57) 

(2.58) 
(2.59) 

(2.60) 
(2.61) 



For T < T c , the recurrence coefficients are asymptotically quasi-periodic, sec DKMVZ]. More precisely, 



lim 



[ln-l{un + ip)} = 0, 



where 



cj = U (T) = 1 - 1 



lim 

n,JV^oo; ^ 



[/3„ - /3( w n + p)] =0, (T < T c ), (2.62) 



R(z) = (z — a 1 )(z - bf)(z - a 2 )(z - b 2 ), 



K = 



dz 



t 1 ^W) 



(2.63) 



7 (x) = -f(x;T), (3{ x) = f 3(x; T) are explicit analytic even periodic functions of period 1 in x, and <p is an 
explicit phase, see ]BDE| ]. The extrema of 7(2;) and j3{x) are expressed in terms of the end-points of the 
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cuts, 



min 7(0;) = 

X 

min [3{x) = 



ai - (a 2 - bi) 



b 2 + a 1 - (a 2 - bi) 



max 7(2;) 

X 

max/3(a;) = 



b 2 - oi + (02 - h) 
4 

&2 + ai + (a2 — bi) 



Using ( |2.4l[ ) , we obtain that as AT — > , 

X = -^ + 0(AT), 



Tr(l-e) 



b 2 ^W) 

hence 

w = e + O(AT) , AT - 
As concerns the extrema of 7(2;) and P(x), they behave as 



2si 



O(AT), 



min, max 7(2;) = 1 ± 

x x 2s\ 



IATI 



1/2 



1 + c 2 



AT + (9(|AT| 3 / 2 ); 



min, max/3(i) = ± — 

xx S\ 



ATI 



1/2 



-^3- AT + 0(|AT| 3 / 2 ), AT^0~ 
24 



(2.64) 

(2.65) 
(2.66) 

(2.67) 
(2.68) 



2.1. Double Scaling Limit for Recurrence Coefficients. We considered above the case when we took 
first the limit n, N — > 00, ^ — ► ^- , and then the limit T — > T c . Here we will consider the double scaling 
limit, when n, N — > 00, 9> — ► 1, with an appropriate scaling of n — AT. We start with the following ansatz, 
which reproduces the quasiperiodic behavior of the recurrence coefficients: 
n 



N 

-,2 



l + N~ 2/3 t, (2.69) 

1 + N~ 1/3 u(t) cos2mre 
+ N~ 2/3 (v (t) + vi (t) cos 2mre + v 2 (t) cos 4mre) 

+ N^ 1 (wo(t) + wi(t) cos 2nire + w 2 {t) cos4ri7re + w s (t) cos6n7re + w±{t) sin4n7re) , (2.70) 
/3 n = + N~ 1/3 u(t) cos (2n + 1)tt£ 

+ TV" 2 / 3 (u (t) cos(2n + l)7re + £ 2 (t) cos (4n + 2)?re) 

+ iV _1 (w (i) + u>i(f) cos (2n + l)vre + iy 2 (i) cos (4n + 2)vre + w 3 (i) cos (671 + 3)?re 
+w 4 (i)sin(4n + 2)7re), (2.71) 

where u(t), Vo(t), . . . ,W4(t) are unknown functions. We substitute the ansatz into string equations ( |2.9| ), 
(2.10) and equate terms of the same order. 

Order iV -1 / 3 . Our ansatz is automatically satisfied at this order. 



Order N- 2 / 3 . We obtain from (EO), (535) that 



1 



-2ci 



-2ci c 2 + 1 

1 -1 
v -l 1 

cf + c| -2cic 2 

-2ciC 2 C? + Co 



By solving these equations we obtain that 



r 2 

c l 2 

4sf 



1 



4s? 



-*T C 





^C 2 


















T 


V-C 3 



T c i A 



and 



Vi - Vi = -u 



v 2 



vo 



4s 2 '' 



Cl 

'4s 2 '' 



r 2 



Is 2 



Cl 

: 2sf 



(2.72) 
(2.73) 
(2.74) 

(2.75) 
(2.76) 
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Order iV -1 . We obtain from (|]|, ( |2.10| ) and fl2.75|) , ( |2.76| ) that 



c 2 + 1 -2ci \ A« \ _uu' ( 



-2ci cf + l! \wq 



1 -l\ /V 



-ci 



-ci 



ci 



8cf 



-1\ T c iu /2cf + 3cf-l 



-1 1 y \vjij 2s 2 \c 2 ) 2s\ \-2c\-4-l 

4 + c 2 . -2cic 2 \ / w 2 \ _ Wi / 1 \ / -2c 2 c 2 

-2c lC2 cj + 4) \ti 2 ) ~ 2 v -c 3 J + 4s 2 V4cf - 3c? + ci 

4 + c 2 -2 Cl c 3 \ ( w 3 \ u^c| f -c 2 + 2 

8s 2 



-C-2 



-4c? 



1 I 4c 2 ' ' 



-2cic 3 c 2 + 4) \w 3 
4 + c 2 -2cic 2 \ / U>4 



-2cic 3 + 1J' 

_ ~2s 2 

-2cic 2 c\ + 4) \wi) 8 ^(4cf - 3c? - 1)J Sl 



(2.77) 
(2.78) 
(2.79) 
(2.80) 
(2.81) 



(we did symbolic calculations with MAPLE). Consider equation ( 2.7q ). The matrix on the left in this 
equation is degenerate, hence we have the compatibility condition, 



2s\u" — u 3 H — % tu, 



(2.82) 



which is the Painleve II equation. When e = 1/2 it reduces to 2u" = u 3 + tu. The function u(t) behaves as 



u(t) ~ — \J —T c t , u(t) ~ Ai(z-ci), k 

t — > — oo Si t — >+oo 



Tc_ 

2s\ 



1/3 



(2.83) 



Here and in what follows we use the following notations: f(x) ~ 5(2;) as x — > a means that lim^^a = 1, 



and f{x) « g(x) as x — > a means that lim :E ^ a [/(a;) — .9(2;)] = 0. 



2.2. Scaled Differential Equations at the Critical Point. Equations (p.7|), (2.8) can be used to derive 
a closed system of differential equations, 



1l — ( ^ n 

N dx \1pn-l 



TV'jx) 



j 2 A n (x) ~fnB n (x) 



\ - 7n B n _!(x) ^^+ 1 lA n {x))\i>n-l 

where 

A„(a;) = x - 4ci + (3 n + (3 n -i, 

B n (x) =x 2 + x(f3 n - 4c x ) + /5g - 4 Cl /3„ + 2c 2 + 7 2 + 7 2 +1 . 
To derive a scaled system at the critical point x = 2c\ we set 

x = 2 Cl +yN~ 1/3 . 

Then 

+ jIMx) = 2 Cl (l - T 2 ) + 7 2 (/3„ + /?„_,) + (7* - 1)^^/3 + ^-2/3 + 
B„(ar) = y 2 + /3„y + 7 2 + 7 2 +1 - 2 - 2 Cl /3„ + ft. 



(2.84) 

(2.85) 
(2.86) 

(2.87) 

(2.88) 
(2.89) 



Substituting ansatz (2.69)-(2.71) we obtain that 
T c V'{x) 



+ jlA n (x) = ^ Cl y2 + Sl^ + S^^ + y UCOS 2nT re -s 1 u'sxn2nne \ + 0(A r ~ 1 ), (2.90) 

T t \ 
y 2 + — + + yu cos (2n + l)7re + silt' sin (2n + l)7rej + C^iV" 1 ) , 

f u 2 T t \ 
7„B n _i = A^-2/3 + _ + + yu cos (2n - l)7re + silt' sin (2n- l)7rej + C^iV" 1 ) . 
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Thus, system (2.84) reduces to 



T, 



d ( ipn \ _ (au(y) a 12 (y)\ ( ip. 



dy \1pn-lJ \ a 2l{y) a22(y)J VV'n-l 



where up to 0(N 1 / 3 ) 



aniv) 


= -ex (y 2 


U 2 

+ y 


0x2(2/) 


= ^ + T 


T c t 
+ 2s{ 


021(2/) 


= -(y 2 + 


u 2 

Y + 


022(2/) 


= ci U 2 + 


u 2 

Y + 



- yu cos 2ri7re + SiU sin 2mre, 
(2n + l)ne + s%u' sin(2n + l)7re, 

T c t\ 

+ —^2 ) ~ yu cos(2?i — l)7re — s\v! sin(2n — l)7re, 



When e = 1/2, this simplifies to 



oix(y) = -022(2/) = -(-i) n y« 5 

,2 



012(2/) = 2/ 2 + ^ + (-iyv, 



021(2/) 



^ 2 +^)+(-i)V 



Under the substitution 



r„(//) : c..s ( n+ 7re/(y) - sin + i ) tt< //( //). 
r„ .j (,7) : cos i n- -) nef(y) -sin(n--) tt< //(//). 



system (|2.9ip reduces, up to 0(iV 1 / 3 ), to 



Si dy \g 



SlU 

y 2 + l 4 + B 



y 2 + u 4 + %)+y"\ ( f 



yu 



the differential ^-equation for Painleve 11 equation ( 2.82) ) 



(2-91) 

(2.92) 
(2.93) 
(2.94) 
(2.95) 

(2.96) 
(2.97) 

(2.98) 

(2.99) 
(2.100) 



(2.101) 



2.3. Universal Kernel. To eliminate dependence on e consider new variables t, u and y such that 



2s\ 



4\ 1/3 



t, u 



4Zc 

Sl 



1/3 



u, y 



47; 

Sl 



1/3 



Then equations (2.82) and ( [2.101 ) reduce to 



and 



d_n 

dy \9 



u" = iu + 2u 3 , (') = --■, 
at 



2v! (Ay 2 + 2u 2 + 1) + Ayu\ // 

(Ay 2 + 2u 2 + i) + Ayu -2u' ) \g 



(2.102) 



(2.103) 



(2.104) 
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Equations ( [2.102 ) give the scaling as 



N 



1 + N~ 2 / 3 



2^ 1/3 



f AT \ ' 

7j 2 = 1 + N~^ 3 f -j^j u cos2nire + 0(N~ 2 / 3 ), 
f AT \ 1 ^ 3 

P n = Q + N- 1/3 ( — -j u cos(2n + l)ire + 0(N- 2/3 ), 



r = 2 Cl +N- 1/3 1 — 

Si 



1/3 



y- 



The Dyson integral kernel for the double scaling limit correlation functions is then: 



K{yi,m) 



/(giMgg) - g(yi)f(m 
yi - V2 



3. Nonlinear Hierarchy 
3.1. Basic Ansatz. For m = 1, 2, . . ., we consider the model critical density 

^2 



p(x) = ^ r (x-2c 1 ) 2m VA~7 2 



where 



2n 



(x-2c 1 ) 2m VA^x 2 dx. 



The corresponding polynomial V(x) is such that 



(.T-2c 1 ) 2m V / 4- x 2 



where Pol[/(x)] means a polynomial part of a function f(x) at infinity. In particular, 



m 



1: V'(x) = — [x 3 - Acxx 2 + 2c 2 x + 8a] , T c = 1 + Ac\ ; 

J C 

2 : V'(x) = — [x 5 - 8cia; 4 + (-2 + 24c 2 )a; 3 - 16cic 2 a; 2 + (-2 + 16cJ 

J C 

+ 16(cj + 4c 3 )] , T c = 2 + 24c 2 + 16c? , 



Be 2 ) a; 



(2.105) 
(2.106) 
(2.107) 
(2.108) 

(2.109) 



(3.1) 
(3.2) 

(3.3) 

(3.4) 

(3.5) 



and so o n. In fact, our considerations will be very gene ral and (3.1) is only an example. They hold for any 
density (1.10) which satisfies regularity conditions ( |1.18| ), ( JL.19 ) everywhere except one point c lying strictly 
inside one of the cuts, and such that as z — > c, h(z) ~ C(x — c) 2m , C ^ 0. 
In the double scaling limit we define variables K, t and y as 



x = j\r~ 1/(2m+1) , 



N 



= l + K 2m Sl t, 



x = 2ci + 2Ky . 



(3.6) 



Our ansatz for the orthogonal polynomials is the following: 

ip(n,x) — cos(n + l/2)nef(t, y) — sin (n + l/2)ircg(t,y) 

+ K [cos (n + 1 /2)7re / x (t, y) - sin (n + 1 /2)?re ^ (t, y) 

+ cos3(n+l/2)7re/(t,y)-sin3(n + l/2)7re5(£,t/) +0(K 2 ), 

[cf. (|2~99| )1, and for the recurrence coefficients, 

7„ = 1 + Ku(t) cos 2n7re + 0(if 2 ) , f3 n = 2Ku(t) cos (2n + l)vre + 0(iT 2 ) 



(3.7) 
(3.8) 
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[cf. fl2.7CD , ( 2.71 )1- See also the work | PeS | where an intimately related ansatz for the recurrence coefficients 
was suggested in the case of a symmetric potential V(x) in the circular ensemble. We substitute ansatz (3.7), 
(3.8) into the 3-terms recursion relation, 

(3.9) 
(3.10) 

(3.11) 



xip(n, x) = j n +iip(n + 1, x) + /3 n tp(n, x) + j n t()(n - 1, x) 
and in the first order in K we obtain two systems of equations, 



(f(t,y) 
\g(t,y) 



L 



f(t,y) 
g(t,y) 



L 



y + u(t) 
-y + u(t) 



(at frequency 1) and 



f(t,y)\ = cMfy (f(t,y) 
g(t,y)J 4s? \g(t,y) 



(at frequency 3). 

3.2. Differential System. We would like to derive a differential equation in y, 



9y 

We are looking for D(t, y) in the form 

D(t,y) = 



f&y) 
g{t,y) 



(f(t,y) 
\g(t,y) 



(3.12) 



-A(t,y) yB(t,y) + C(t,y) 

yB(t,y)-C(t,y) A(t,y) 



(3.13) 
(3.14) 



[cf. ( 2.101 )1, where A, B and C are even polynomials in y of the following degrees: 

degA = 2m-2, degB = 2m-2, degC = 2m. 

We will assume that C is a monic polynomial, so that C = y 2m + . . .. The general case can be reduced to 
this one by the change of variables, t = nt, y = , u(t) = , which preserves the structure of the operator 
L in (3.1C). The consistency condition of equations ( |3.10| ) and (3.12), 

1' 



[D, L] = d y L - d t D = 



-1 



dtD, 



(3.15) 



implies that 

d t B = 2A, d t C = l + 2uA, d t A = ~2y 2 B + 2uC . (3.16) 
Example: m = 1. According to ( |3.14 ), A = a(t), B — b(t), C = y 2 + c(t). From the last equation in 



(3.16) we obtain that b = u and then that 

u' 

a = 2 

where to is a free constant, and 



c = t+ — +t a , 



u 

T 



u 3 + 2(t + t Q )u, 



(3.17) 



(3.18) 



the Painleve II equation. By changing t + to to t we can reduce it to to = 0. 

We would like to construct solutions to ( 3.16 ) for m > 1. To that end, define recursively functions 
A m (t,y), B m (t,y), C m (t,y) by the equations 

C m+1 = y 2 C m + f m (u) , (3.19) 

B m+1 =y 2 B m + R m {u), (3.20) 



1 



A m +i = y A m + - d t R m (u) , 



where R m (u), / m (it) solve the recursive equations 



R m +i(u) = uf m (u) - - d tt R m (u), 
d t fm{u) = ud t R m (u), / m (0) = 0, 



(3.21) 

(3.22) 
(3.23) 
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with the initial data 



,40 = 50 = 0, C = l, Ro(u)=u, f (u) = — . 



We solve recursively ( |3.22[ )-(|3.24D as 

Ri(u) = «", h(u) = |« 4 - \ uu" + \u' 2 

R 2{ u)^~\u^u"~\uu' 2 + ^\ 



f 2 ( u ) = — u 6 - - u 3 u" ~ — u 2 u' 2 + —uu^ - —u'u" 1 + —u" 2 



■5 

16 



1 

16' 



16 



1 

32 



n i \ 5 7 35 4 ,, 35 o ,2 , 7 2 C4~i , ' / /// , 21 „2 , 35 ,2 1 ir\ 

R 3 (u) = — u' u u u u H u u ( ' + - uuu H uu H u u u w , 

iK ' 16 32 16 32 8 32 32 64 



and so on, 



In addition, 



R m (u) 



(2m)! o m+1 



2 2m (m!) 2 



(~l) r 

22m 



,(2m) 



A x = - £1 = u, Ci = y 2 + i u 2 ; 

A 2 = - u'y 2 + - uV - - u"', B 2 =uy 2 + -u 3 -- u" 
2 J 4 8' y 2 4 



(3.24) 

(3.25) 
(3.26) 
(3.27) 
(3.28) 

(3.29) 

(3.30) 
(3.31) 
(3.32) 



and so on. It is easy to check that the functions A m (t, y), B rn (t,y), C m (t,y) defined by ( |3.19 )-(3.24) solve 
the equations 

d t B m = 2A m , d t C m = 2uA m , d t A m = -2y 2 B m + 2uC m -2R m (u). (3.33) 
Indeed, by fl3.24| ) it holds for m = 0. Assume that it holds for some m. Then by (|1|)-(U|) and (|3|), 

(3.34) 
(3.35) 



dtC„ 



m+1 



y 2 d t C m + d t f m (u) = 2y 2 uA rn + ud t R m {u) = 2uA 
d t B m+1 = y 2 d t B m + d t R m (u) = 2y 2 A rn + d t R m (u) = 2A m+1 , 

d t A m+1 = y 2 d t A m + i d tt R m {u) = y 2 [-2y 2 B m + 2uC m - 2R m (u)] + - d tt R m (u) 

= -2y 2 [y 2 B m + R m (u)] + 2u [y 2 C m + f m (u)] - 2 (uf m (u) - j d tt R m {u^j 

= -2y 2 B m+1 + 2uC m+ i - 2R m+1 (u) , (3.36) 
which proves ( [3.33 ) for m+1 and hence for all m = 0, 1, 2, . . .. Comparing (3.16) with (3.33) we obtain that 

A = A m , B = B m , C = t + C m (3.37) 
solve equation (|3.16 ), provided u(t) is a solution of the equation 

R m (u)+tu = 0. (3.38) 

The sequence of equations ( |3.3£ ) for m = 1, 2, . . . forms a hierarchy of ordinary differential equations which 
is known as the Painleve II hierarchy [Kit] (see also [Moc], [PeS|). We can formulate now the following 
result. 

Theorem 3.1. Define 

n (+ ,,\_( -A m (t,y) yB m {t,y) + C m (t, y y 

- {y Bm (t,y) - C(t,y) +A m (t,y) 

Then if u(t) is a solution of equation ( 3. 3t\ ), then the matrix 



(3.39) 



D(t,y) 



t 
-t 



D m (t, y) 



(3.40) 
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is a solution to (3. It). More generally, if ti 
equation 



, t m are arbitrary constants and u(t) is a solution of the 
Y,tkRk{u)+tu = 0, (3.41) 



fe=i 



then the matrix 



D(t,y) 



t 
-t 



(3.42) 



fc=i 



is a solution to ( 3.11 ). 



Remark: It can be shown that ( 3.42 ) is a general solution to equation ( 3.15 ). 

The meaning of the constants ti, . . . ,t m in ( 3.42| ) is the following. Observe that the differential equation 
in y, (3.12) describes the double scaling limit for a critical polynomial of degeneracy 2m. In this case the 
space of transversal fluctuations to the manifold of critical polynomials has dimension m. The variables 
ti, . . . , t m serve as coordinates in the space of transversal fluctuations, and ( |3.42| ) gives the matrix describing 
the double scaling limit of the recurrence coefficients in the direction r = (ti, .. . , t m ). 

4. Conclusion 

In this paper we considered critical polynomials which violate the regularity conditions at exactly one 
point, inside the support of the equilibrium measure. It is characterized by the degree 2m of degeneracy of 
the equilibrium density at the critical point. Our main results are the following: 

• When m = 1, the infinite volume free energy exhibits the phase transition of the third order. This 
extends the result of [GW] to nonsymmetric critical polynomials. 

• When m — 1 , the double scaling limit of the recurrence coefficients is described, under a proper 
substitution, by the Hastings-McLeod solution to the Painleve I I diff erential equation. Befor e thi s 
result was known only for symmetric critical polynomials [ DSS ], |PcS| (for rigorous results see |BI2|, 



• For m > 1, we derive a hierarchy of ordinary differential equations describing the double scaling 
limit of the recurrence coefficients. 

5. Appendix. One Useful Identity 



Let V(z; T) = -^fc^ , where T > is the temperature, and 



fi N (dM;T) = ZjvCTj^exp 



N 



Tr V(M) dM, Z N (T) 



exp 



TV 
T 



Tr V(M) dM . (5.1) 



Then p(x) and u>(z) depend on T. The following identity is useful in many questions. 

Proposition. Assume that the number of cuts does not change in a neighborhood of a given T > 0. Then 



d_ 
~dT 

where the numbers bj < Xj < aj + i, j = 1, . 



RV 2 (z) ' 
q — 1, solve the equations, 



ak+1 ULlix-xj) 



dx = 0, k = 1, 



R}/ i {x) 

The neighborhood can be one-sided, then the derivative in T is also one-sided. 
Proof. Equation (|Tl|) gives that 

V'{z) h{z)R 1 ' 2 {z) 



Since V{z) does not depend on T, 



Tlo{z) 



d d h(z)R^(z) 
dT [Tuj{z)] = -dT 2 



(5.2) 



(5.3) 



(5.4) 



(5.5) 
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The function on the right can be written as 

d Th(z)R 1 l 2 {z) P(z) 

dT 2 R l / 2 (z) ' 

where P{z) is a polynomial with real coefficients. Since 



we obtain that 



^;[Tuj(z)} = - z +0(z- 2 ), 



which shows that P(z) = zJ- 1 + .... By fll.l5| ), 

n +i h{x)B}i 2 {x) 

L 2 



dx = , j = 1, . . . ,q - 1 . 



By differentiting with respect to T we obtain that 

P(x) 



dx = , j = 1 , . . . , g — 1 . 



This is possible only if P{x) has a zero in each interval [bj, Oj+i]. Thus, (5.3) is proven. 
As a corollary, from (|T^) we get that 

m&z-xj) 



dT 



h(z)R 1/2 (z) 



This implies that 



— In 

dT 



h(z)R 1/2 (z) 



R}/ 2 (z) 



h{z)R{z) 

Comparing the residue of the both sides at z — a k ,b k we obtain that 



dak _ 

dT h{a k ){a k -b k )Y\ j . :j _ tk [(ak- a ] ){a k -b j )] ' 

dbk = ^UjZliak - xj) 

dT h(b k )(b k - a k ) Ur.^kKh ~ a,)(b k - bj)} ' 



1 < k < q. 



(5.6) 

(5.7) 
(5.8) 

(5.9) 

(5.10) 



(5.11) 



(5.12) 



(5.13) 



[BDE] 
[BDJ] 

[BI1] 

[BI2] 

[BIPZ] 
[BIZ] 

[BKa] 
[BPS] 

[BZ] 

[CDM] 

[DGZ] 
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